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Abstract 

Using perturbed Bogoliubov equations, we study the linear response to a weak orbital mag- 
netic field of the heat conductivity of the normal-superfluid interface of a polarized Fermi gas at 
sufficiently low temperature. We consider the various scattering regions of the BCS regime and 
analytically obtain the transmission coefficients and the heat conductivity across the interface in 
an arbitrary weak orbital field. For a definite choice of the field, we consider various values of the 
scattering length in the BCS range and numerically obtain the allowed values of the average and 
species-imbalance chemical potentials. Thus, taking Andreev refiection into account, we describe 
how the heat conductivity is affected by the field and the species imbalance. In particular, we show 
that the additional heat conductivity due to the orbital field increases with the species imbalance, 
which is more noticeable at higher temperatures. Our results indicate how the heat conductivity 
may be controlled, which is relevant to sensitive magnetic field sensors/regulators at the interface. 
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I. INTRODUCTION 



Recent experiments on ultracold Fermi gases [l-^, such as and ^Li atoms, have 
shifted our perspective on superfluidity and confirmed earher ideas about fermionic paring 



6|. The extent of possible experimental control in these gases has led to a large body of 
experimental work on fermionic pairing phenomena. Most notable is the ability to apply 
a controlled magnetic field to tune the strength of the attractions 2|, |3|. In balanced mix- 
tures of fermions, this tunability is exploited to study the crossover from a Bose-Einstein 
condensate (BEC) of molecules to a Bardeen-Cooper-Schrieffer (BCS) superfluid The 
BCS regime involves a small negative scattering length wherein fermions form weakly bound 
Cooper pairs, while in the BEC regime the scattering length is positive and particles form 
tightly bound molecular pairs. However, in the case of imbalanced mixtures of spin-polarized 



Fermi gases 
experimental 



7 



le study of the phase diagram is more interesting 0, [lo|. In this case, many 
ll|-|l4| and theoretical 15|-|26| studies have been reported. In the experi- 



ments, a controllable number of atoms prepared in one spin state can transfer to the opposite 
spin state by applying a radio-frequency pulse. In these polarized systems, all atoms in the 
less populated spin state are found to pair up with atoms in the opposite spin state, consti- 
tuting a molecular condensate that coexists with the remaining Fermi sea of the unpaired 
atoms. As the repulsive interaction between these atoms and the molecules increases, a 

n 

phase separation occurs between the superfluid (SF) and the normal (N) phase [7[. Such a 
phase-separation scenario had been proposed by Clogston 271] and Chandrasekhar 28|] long 
ago, who predicted the occurrence of a flrst-order transition from the N to the SF state. Since 
the existence of a N-SF interface may settle some important questions regarding polarized 



Fermi mixtures, serious research has been focused on this subject 
it has been established 



29 
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33| . In particular. 



31| that a temperature difference between the two phases can 



appear as a consequence of the blockage of energy transfer across the N-SF interface. This 
blockage is due to a SF gap, which causes low-energy normal particles to be reflected from 
the interface. By studying particle scattering off the interface, the heat conductivity has 
been also calculated. We may mention in passing that the heat conductivity has been cal- 
culated in other polarized Fermi systems, such as the spin-polarized superfluid Fermi liquid 



32], where a phase separation does not occur. 

In this article we examine the N-SF interface of a polarized Fermi gas in the presence of a 



2 



wea. 
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i orbital magnetic field, which is produced fictitiously, for example, by rotating the gas 



35| . Under rotation, the neutral system behaves like a charged system in a magnetic 



field. Another way to produce fictitious magnetic fields is via vortex nucleation [36|. Also, 
when a neutral system moves in a suitably designed laser field, it can mimic the dynamics of a 
charged particle in a magnetic field [34] . We use perturbation theory to solve the Bogoiubov 
equations and obtain the wave functions. We consider the various scattering regions of the 
BCS regime associated with different energies of the incident particle/hole and analytically 
obtain the transmission coefficients and the heat conductivity across the interface. For 
definite choices of the field, considering various values of the scattering length in the BCS 
range and taking Andreev refiection into account, we describe how the heat conductivity is 
affected by the orbital field and the species population imbalance. In particular, we show 
that the additional heat conductivity due to the orbital field increases with the species 
imbalance at higher temperatures. Our results indicate how the heat conductivity may be 
controlled, which is relevant to sensitive magnetic field sensors/regulators at the interface. 



II. PERTURBED BOGOLIUBOV EQUATION 

We consider a Fermi gas consisting of two imbalanced species a and b, each with mass 
m, effective charge e and chemical potential /ij {i = a,b). (The external trapping potential 
is assumed to be zero.) The Hamiltonian operator in the presence of the orbital magnetic 



field is = 1) |37| 



Ha,b = -7^{V-teAY-fiTh (1) 
2m 

where = {fia + /^fe)/2 and h = {fia — /ib)/2 are the average and the species-imbalance 
chemical potentials, respectively. Treating the vector potential as a perturbation, we may 



use the first order perturbation theory to solve the Bogoliubov equations 38|. Denoting the 



particle- like and hole- like solutions by Ui{r) and Vi{r), respectively, we have 



U4 = u 



(°)+««, y^ = ^ + ^^ a = AW+A(i) (2) 



where uf\r) and vf\r) are the solutions to the unperturbed (A = 0) Bogoliubov equations 
and A(r) is the pair potential. We use the approximation that A*^") is independent of r and 
by suitable choice of A (see later) the contribution of its component perpendicular to the 
SF boundary can be made negligible. Thus, and using the London gauge (V ■ A = 0), one 



finds that A^^^ ^ [39]. The Bogohubov equations then read: 

( + E + fi + h]u^^^- Avi''^ = -A ■ Vwi°) 

\2m r~ J a " m 

(-^V^ + E-^^ + h) - = . Vvf\ (3) 

V 2m / m 

The second set of Bogoliubov equations is obtained by simply interchanging a and h. It 
is clear from these equations that the incoming a-particles are coupled to 6-holes and vice 
versa. To proceed, let us take the N-SF interface to be in the x = plane. Introducing the 
superscript s (n) for the solutions in the SF (N) phase, we expand u\^^^\vl'^^^^ in terms of 
the unperturbed eigenf unctions. The latter are (jy^g-^ir) = exp[i(fc|| ■ r ± k(q)x)] for the N 
phase, and (j)k^''\r) = exp[i(fc|| ■ri/c^'^^x)] for the SF phase, where q = p,h refers to particle, 
hole and fey {k^'^\ k(^g)) denotes the component of the wave vector k parallel (perpendicular) 
to the interface. Notice that q appears as a subscript (superscript) for the N (SF) phase 
throughout our notation. From the unperturbed Bogoliubov equations we have the usual 
relations 

k^h) = ,/^~~^^^^^{ETli}, ^^"''^ = ^kl) - 2me^ (4) 



where ^± = E + h ± ^ {E + hy — A"^. Thus, for the N phase, we write 
wlfeV) = f^fc(p)0fc(p)(^) + perturbation terms (PT) 

cr=± 

= E + perturbation terms (PT). 

cr=± 

As for the SF phase, 

q,a q,cr 



u 



ak 



E C^'Vr (r-) vlT ^(r) = E CVr (^) (6) 



k',q,cr k',q,(T 



where V^"-"' = BU^^^' and V^^'"' = B-^U^^""' with B = ^_/A. Substituting (E]) into (ED, we 
can use the orthonormality of 0's to extract the coefficients a^}^) and b'^^^) . We obtain, after 



some calculations 



where 

= E + /i±[/.-^(A;J + fc(^)^)]. (8) 

In our calculations we have naturally supposed that for cr 7^ a', 

/ d'x [C^'^^^VC'^") - ^'("^0?)*] ■ A = 0. (9) 

(The integrand, then, does not pertain to a single entity (particle/hole) and is thus irrelevant 
as far as the transmission of particles/holes is concerned.) The amplitudes U and V are to 
be determined by matching the wave functions and their derivatives at x = 0, of course. 
This is deferred to the next section. 



III. TRANSMISSION COEFFICIENTS 

Consider an incoming a-particle from the N side with energy E > C,{p) — l^a-, where 
^(p) = ^p^/2m. Because there is no correction to the energy to lowest order in A 3^, in 
view of (jl]), there are different scattering regimes involved in the E — C,{p) plane. These 
regimes are obtained by the crossover of either of the k^^-l^"^ or fey from real to imaginary 
values (or vice versa), signifying a change in the scattering mechanism as described below 



3G, 



3l|: 



i) For E < A — h, the incoming particle has insufficient energy to excite the SF side. In this 
case we have total reflection. 

ii) For E > A — h and C,{p) < the incoming particle is reflected. 

iii) For E > A — h and < ,^(p) < only quasi-particles are excited in the SF. 



iY)FoT^+<^^p)<2{E + h 
but Andreev reflection 



particle-like and hole-like excitations both occur in the SF side 



40l442l| is forbidden. 

v) For 2{E + h) < C,{p) < E + fia, we have particle-like and hole-like excitations as well as 
normal and Andreev reflections. 

Because in regions i and ii the particle has insufficient energy to excite the SF side, the 
transmission coefficients in these regions vanish. Moreover, our focus is on energies slightly 
above the transmission threshold {E ^ A — h), because we are considering low temperatures 
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(see later). Therefore, we concentrate on regions iv and v. We proceed to obtain U and V 
by matching the wave functions and their derivatives. For region iv we have 



where 



1 + I]fc' ^kk' 



1 + Yjk' ^kk' 

( 









k(p) 



k(p) 



-iP) 




PT 



k(h) 



PT 



X, 



kk' 



M, 



-(h) 



XL' = M 



kk' 



-ip) 
kk' 



A2 



jiP) j(p) _ ^2 

fP + AB 



kk' 



(10) 



(11) 



Similar equations are obtained for region v. These equations yield the scattering amplitudes 
U'^^''^ for general A. Denoting the x-component of the current density by j^, the transmission 
coefficient is given by = jj/j^, where the superscripts T and I refer to the transmitted 
and incident quasi-particle current densities, respectively. The general form of j (for channel 
a) is 



jc.{r) = - — [<(r)Vna(r) - Ma(r)V<(r) - vl{r)Vvb{r) + Vb{r)Vvl{r)] 

--[<(r)«,(r) + vl{r)v,{r)]A. (12) 
m 

Using 05]) and we obtain, to first oder in the perturbation, 

W = + Uvikk'^'^ - k^'^) + v'kl'ik'^'^ - k^'') + Vi^lik"^ + k^'^) + r/S(fc'(^) + A:(^))] 



"(p) k' 



(13) 



where W^^^ is the transmission coefficient in the absence of the vector potential and 



jj+iP>^-ih) 



u+ip)u-w* 



(1) _ ^k ^k Y \ r,(2) _ ^k Y' \ 

Vkk' - 77X2 [^kk' - ^kk'), Vkk' - 77+2 y^kk' - y^kk') 



k{p] 



'^Hp) 



„(3) _ ^ 
'Ikk' 



-{h? 



B '^Ykk' - Xkk'), Vkk' 



(4) _ f^fc 



<P? 



k{p) 



(Vkk' - B XL')- 



(14) 



Hp) 



In the computations (which are performed analytically throughout) we use 



l7^s=—JdKldk[^, (15) 

k' 



(27r)3 Att^J "7 

as we are considering the a channel. Note that only energies slightly above the transmission 
threshold contribute to the first integral, which affects the region of integration. We also 
need to make a definite choice for A. Equation (flQj) then yields very complicated analytical 



expressions for the scattering amplitudes that reduce to the results of 311] in the absence of 
the vector potential. Substituting these in (I7j) gives the coefficients a'^^^) and b'^]i^) . 

We take A = — + , ^3 7" + , .^ k, where c and xn 7^ are constants. (This 
particular choice renders analytical computation possible.) Equation ( IT3|) leads to 

W, = wi''^2ec,I^^^Z±±^^(A-k)^E^AT-h (16) 

V X XT VX^ - A2 

where x = ^ip)~^ ^+ (^-) refers to region iv (v). coincide with the results of jsil 
in the absence of the vector potential (c = 0), which are also proportional to \/ E — A + h. 
We remark that only a tiny portion of region iv (corresponding to C,{p) slightly above 2 A 
and, of course, E slightly above A — h) contributes to W^. Similarly for holes, we obtain 
results for the region v with C,{h) = ^fh)!'^'^- This is because from our calculations, holes 
incident on the interface can only transmit in region v. 

Analogous procedure applied to incoming 6-particles and a-holes (the (3 channel) yields 
identical results to the above but with h ^ —h. 

IV. HEAT CONDUCTIVITY 

For each channel, the heat conductivity is given by 

^ = E / d^M I dEEmWiE,^^,)) (17) 

q=p,h 

where f{E) is Fermi-Dirac distribution, which, for sufficiently low temperatures T << A 
(^B = 1), reduces to e~^^~'^'>^'^ (up to a proportionality constant). This implies that the 
relevant range of energies is only slightly above the threshold E = A — h, a.s mentioned 
before. 

Using (dZD we find k = + k^^) = (k(^°) + k^^) + (/t^ + k^), where is the heat 
conductivity in the absence of the external field and k*-^-* is the additional term due to the 
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FIG. 1: (Color online) Heat conductivity versus temperature across the interface in the BCS regime 
for c = 0.45. The dotted line represents the heat conductivity in the absence of the vector potential. 

vector potential in the SF side. The effect of the vector potential in the N side on k^-^^ is of 
second order and, therefore, negligible. Its effect on k^^^ is to add a small term proportional 
to K^^^ itself. This does not change the temperature dependence of k^^^ and only causes a 
small shift in the k^^^ — T curve (figure 3 of reference 30|), which we neglect. Also, in our 
calculations we have neglected the energy carried by the /3-channel because it is quite lower 
than that carried by the a-channel. Our final expression for is 



K 



(1) 



12v/2ecA(A - hf/^m^/\V2A - J fi + A + 



1 -{A-h)/T 

) 



(18) 



By fixing the values of h and /i using the gap and number equations [31|, |43j, we plot 
in FIG. [H k/kn (/tN = 2m/iT/7r^, heat conductivity in the normal phase) as a function 
of T/Tf in the BCS regime (A = 0.25Tf and (/cpa)"^ = -0.44). The effect of the vector 
potential is simply a shift of the heat conductivity curve. We have plotted the figure in 
the range O.OSOTp < T < 0.045Tf, where the shift is more pronounced. Also in FIG. |2]we 
plot K^^y versus T/Tp, for various values of c. As seen, at fixed T, the larger the value 
of c, the larger is the heat conductivity. However, at temperatures roughly below O.OlTp 
the effect of the vector potential is negligible. Moreover, the curves show an increasing heat 
conductivity, with a rate that decreases substantially after around T = O.OSTp. FIG.|3]shows 
the dependence of k*^^)/^n on /i at fixed temperature for various values of T/Tp. To obtain 
this plot, we have first fixed Aa^ (a, the scattering length). Using the gap and number 
equations, we then fix kpa (such that it lies in the BCS range) which fixes the values of /i, 
h and A. As seen, in the specified temperature range, n^^^ does not change significantly 
with fj, at lower temperatures but decreases as the temperature is raised. Notice also the 
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FIG. 2: (Color online) Additional heat conductivity versus temperature across the interface in the 
BCS regime for various values of c. 
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FIG. 3: (Color online) Additional heat conductivity versus average chemical potential across the 
interface in the BCS regime for c = 0.1. 

(almost) collapse of the data points, showing that the graph is roughly independent of the 
temperature as long it is high enough. We also show, in FIG. HI the dependence of k^^^ / on 
h at fixed temperature for various values of T/Tp. Again, k^^^ does not change significantly at 
lower temperatures but now increases (in contrast to the previous figure) as the temperature 
is raised. Similarly, the graph is roughly independent of the temperature for high enough 
values. Since h can be controlled by the species population imbalance, the graph shows how 
the effect of the weak magnetic field on the heat conductivity may be controlled. 
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FIG. 4: (Color online) Additional heat conductivity versus imbalance chemical potential across the 
interface in BCS regime for c = 0.01. 
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